Intrinsic time gravity and the Lichnerowicz-York equation 
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We investigate the effect on the Hamiltonian structure of general relativity of choosing an in- 
trinsic time to fix the time slicing. 3-covariance with momentum constraint is maintained, but the 
Hamiltonian constraint is replaced by a dynamical equation for the trace of the momentum. This 
reveals a very simple structure with a local reduced Hamiltonian. The theory is easily generalized; 
in particular, the square of the Cotton- York tensor density can be added as an extra part of the 
potential while at the same time maintaining the classic 2 + 2 degrees of freedom. Initial data con- 
struction is simple in the extended theory; we get a generalized Lichnerowicz-York equation with 
nice existence and uniqueness properties. Adding standard matter fields is quite straightforward. 
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INTRODUCTION 

The Hamiltonian theory for general relativity was most 
clearly laid out by Arnowitt, Deser, and Misner [l[ more 
than 50 years ago. The phase space consists of a pair 
(fe , tt ij ) , where is a Riemannian 3- metric and tt 13 is 
the conjugate momentum. These cannot be freely chosen 
because they must satisfy the constraints 
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til TT 
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0, 
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where tt = gijir^ is the trace of tt 1 K These conditions are 
known respectively as the Hamiltonian and momentum 
constraints. We must also choose a scalar and a vector 
(N,N l ). These are the 'lapse' and 'shift' . The total 
Hamiltonian, in the compact without boundary case, is 



\tt 2 ) - 2N j V i irV]d 3 x. 



H = J [g- 1 / 2 N(-gR + ir ij 'Ki. 

(2) 

Therefore (N,N l ) are the Lagrange multipliers of the 
constraints. The first of Hamilton's equations gives the 
relationship between tt ij and the time derivative of gij 



dt 



2Ng- 1 l 2 {^ ] ~ ^ 9ij n) + V,Y, + V, N t . (3) 



A major difficulty with the canonical quantization of 
gravity program is this freedom to choose -/vQ. Each 
choice of N gives a different slicing of spacetime, and 
reflects the 4-covariance of the Einstein equations. We 

want to break this covariance, choose a natural time vari- 

i — I 

able, and compute the emergent lapse. James York in |3| 
pointed out that the local volume, ^fg, and tt j ' ^fg are 
canonically conjugate, and suggested that tt / ' ^fg is a nat- 
ural time, and that ^/g is the local energy. This is called 
an 'extrinsic' time because 7r y is essentially the extrinsic 



\gijv) 



K, 



curvature of the slice. We have g~ 1 / 2 {TT 
where is the extrinsic curvature. With York's choice, 
^fg is the associated energy density and the total volume 
of the slice is the Hamiltonian. Two years earlier Charles 
Misner |4j made the opposite choice, albeit in a minisu- 
perspace context. He picked the local volume as his time 
and 7r is then the energy density. This is where he in- 
troduced the 'mixmaster universe'. The great advantage 
of Misner's choice is that a local reduced Hamiltonian 
appears. Recently, two of us [2| presented a theory of 
gravity passing from classical to quantum regimes with 
a paradigm shift from 4-covariance to 3-covariance. Its 
framework revealed the primacy of a local reduced Hamil- 
tonian and intrinsic time proportional to lng 1 / 3 . The 
present article is an expansion and extension of part of 
that work. 



INTRINSIC TIME GRAVITY 

We start by choosing the intrinsic time function as 
lng 1 / 3 . Therefore Sing 1 / 3 , for any variation, is a scalar. 
Further, the canonical conjugate is tt, with coefficient of 
unity. We split the metric into unimodular and determi- 
nant parts via gtj — g^^gij, and split tt 1 '^ into tracefree 
and trace parts via 7f y = g 1 / 3 ^^ — g^ir/i). The sym- 
plectic 1-form becomes 



ttS In g 



1/3 



(4) 



We see immediately that (gij,^ ) and (lng ' ,tt) form 
conjugate pairs, which is clearly the generalization of the 
PiSx 1 + PoSt one gets in the case of a simple particle. 
With the choice of time as In g 1 / 3 (which varies from — oo 



2 



to +00, instead of to 00), we can see that tt is analogous 
to P = -E. 

We now substitute the decomposition of (gij,TT t: >) into 
the Hamiltonian constraint, Eq.([T]), to give 



gR + 9ik9jiTr ij n kl - (3' z tt 



0. 



(5) 



where f3 2 = | for GR. It is a free positive parameter for 
the extended theory which we discuss later. We know 
that the Hamiltonian as the generator of time trans- 
lation, is conjugate to time and should equal the en- 
ergy, and therefore the true Hamiltonian density with 
this choice of intrinsic time should be —tt. We then just 
solve the Hamiltonian constraint, which is equivalent to 
(tt - ff/ v /^2)( 7 r + HjsfW) = 0, to find the reduced 
Hamiltonian 



B^ = -^g lk g^TT kl 



gR; 13 = ±y/p. (6) 



A simple toy model for this process is given by the 
relativistic particle, which satisfies a constraint 
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(7) 



P is the energy, and therefore the physical Hamiltonian 
is 



E = P 







-P = H 



P P 



(8) 



Hamilton's equations give the equations of motion for 
P. We then add the constraint P° = H as a dynamical 
equation to determine P°. 

The Hamiltonian H/0 of Eq.© generates In g 1 / 3 {x,t) 
translations, and gives equations of motion for g~ij and 
7f y with respect to this intrinsic time variable. We need 
to write gR in terms of cjij and lng 1 / 3 . This is quite 
straightforward. We do not have an equation for Ing, 
since it is the time, but we do need a dynamical equation 
for 7T which is given by 



0. 



(9) 



This is a rewriting of the Hamiltonian constraint, but it 
is no longer to be viewed as a constraint; rather, it is 
the evolution equation for tt. This is the fundamental 
equation for intrinsic time gravity; and it can be inter- 
preted respectively as the Hamilton- Jacobi equation and 
Schrodinger equation in the semi-classical and quantum 
regimes [2|. 

We choose the positive root for H when we take the 
square-root. Negative values of tt correspond, in line with 
current observations, to an expanding universe. The re- 
duced Hamiltonian can be constructed before solving for 
tt. Since the system is 3-covariant, the evolution equa- 
tions will propagate the momentum constraint. This is 
the only constraint left since the Hamiltonian constraint 



is gone because we have a unique choice of (intrinsic) 
time. Thus, we are left with a system that has the ex- 
pected 2 + 2 degrees of freedom. 

The symplectic potential of the conjugate pair 



(lng 1 ' 3 ,^) contributes / J(n 



ding 1 
dt 



-)d 3 xSt 



~ III TF aiU dt 1 d 3 x]dt, to the action; and it is thus 

clear that J jf dln g t ' d 3 x contributes to the total Hamil- 
tonian generating ^translations, where t is the ADM 
time parameter. By subtracting the tangential change 
generated by the the momentum constraint, the rate of 
change of the normal component of In g 1 ^ 3 is 



lim 

st^o 



Sing 1 / 3 
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(10) 

The classical evolution of (fl , i_/,7r y ) w.r.t. the ADM time 
variable t can equivalently be obtained from the effective 
Hamiltonian 

Hf = J [f(&/0) + N*Hi] d 3 x 

=1 [fiy/mg^'*" 1 -gR- 2N l ( g -^v^)]d 3 x(ii) 

On the other hand, by contracting Eq.© with g % \ we 
find 



Sing 1 / 3 
dt 



If we choose 



N 



Ntt 2 
Wl + 3 



(12) 



(13) 



it is clear that Eq. tfTU]) agrees with Eq. (|T2"j) . and if one 
uses the lapse as given by Eq. fTBI . it is a straightforward 
exercise to show that the evolution equations arising from 
Hf agree with those from the ADM evolution equations. 
Given any t- foliation of spacetime, we can find a corre- 
sponding reduced Hamiltonian that generates it. This 
means that 'many fingered time' lives on in the intrinsic 
time picture but in a very different form by the way the 
emergent N, as given by Eg . ([13]) . depends on the (freely 
chosen) /. 

We have expressed GR in a very simple form: we have 
only the effective Hamiltonian, Eg . (fTTj) . combined with 
the dynamical equation for the trace of the momentum, 
Eq.©, and the momentum constraint. 

Moving away from vacuum GR, how can we generalize 
this structure, while maintaining the 2 + 2 degrees of 
freedom and the 3-covariance? 

There are two changes we can make. We can multiply 
each of the three terms in the Hamiltonian constraint, 
as given by either Eq.([T]) or Eq.([5|), with arbitrary con- 
stants. This means deforming /3 in the tt 2 term away 
from jr, which is equivalent to changing the constant in 
the DeWitt supermetricQ. We really only need two con- 
stants because we can always rescale one of them to unity. 
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Thus we replace gR by a 2 gR and leave the gikSjl^^ 
alone. More radically, we can replace the 'potential', —R, 
by any scalar function V of the metric, and everything 
still works. The set of equations now reads 

H t = I ($y/9ik9ji* ij * kl + gV + N'H^cPx, (14) 

^ := -^ySikSjl^^ 1 + gV, (15) 
da 1 / 3 2 

/ := - V - 3 V^. (16) 

INITIAL DATA AND EXTENSIONS OF THE 
LICHNEROWICZ-YORK EQUATION 

It is not quite that easy, however; we need to find ex- 
plicit initial data. There is no point in having a nice 
Hamiltonian system with no solutions. We need to find 
a pair (gij , ir^ ) that satisfies both (3tt + H = and 
Vi7r iJ = 0. In the case of GR this pair reduces to 
the standard constraints; and the only really successful 
general way of solving them is the conformal method, 
which was initiated by Lichnerowicz [||. There is a very 
comprehensive account in Q, especially in Chapter VII. 
The technique is to choose free data that consist of a 
base metric g^, a symmetric tensor density tTj TT that 
is both tracefree and divergence-free with respect to g, 
and a scalar p. It is particularly simple if p is a con- 
stant. This guarantees that the extrinsic curvature has 
constant trace, and thus we construct a 'constant mean 
curvature' (CMC) slice. We make a conformal transfor- 
mation g^ = 4> A gij, 7r u = <f> tt 13 . It turns out that Tr* TT 
is conformally invariant; it remains TT with respect to 
gij. Thus we write 7r- TT = ^} TT - We also set n = y/gp. 
We are guaranteed that 7r IJ = (njl T + y/gg^p/S) satis- 
fies the momentum constraint for any conformal factor 
<t>. In particular, we can specialize to (j) 4 — g 1 ! 3 to get 
(g~ij,ir l: >) = (g^, #y T ). We now seek an appropriate in 
order to solve the Hamiltonian constraint. In GR, this 
reduces to solving the Lichnerowicz- York (LY) equation 

8V 2 -Rcj> + T 7rf V 7 ~ \?<f = 0, (17) 

where g is the determinant of gij. This is an extremely 
nice equation because it always has a unique, positive 
solution 0. 

Let us now multiply the R and the ir 2 terms in the 
Hamiltonian constraint by arbitrary positive constants. 
This means replacing Eq.([5]) by 

- a 2 gR + g.L.gji-"- 1 ' 1 ~ [3 2 ir 2 = 0. (18) 

The new LY equation becomes 

8a 2 t 2 (p-a 2 R<f> + g- 1 Ttf T 7Tj TT <p- 7 -/3 2 p 2 <p 5 = 0. (19) 



This equation is just as nice as the original LY equation, 
Ec flTTj) . because it too always has a unique positive so- 
lution. Such a rescaling has been recently discussed in a 
different context in 

However, if we add any other function of the metric to 
Eq. ([Tg)) we destroy the nice properties of the LY equa- 
tion. For example, if we were to add an R 2 term, the LY 
equation would pick up the term (8V 2 — R(f>) 2 , which 
changes the nature of the LY equation completely and 
the existence and uniqueness results no longer hold. 

There is one exception. York, in [l(|, rediscovered a 
conformally covariant tensor which is a function of the 
metric. This is now known as the Cotton- York ten- 
sor(density), i.e., /3 lJ . This transforms exactly as a TT 
tensor does under conformal transformations, so /3j is 
conformally invariant. Therefore we can pick a general- 
ized Hamiltonian as 

- a 2 gR + g ik 9jiTt ij Tt kl + 7 2 /3j^ - (3 2 tt 2 = 0, (20) 

where 7 is another coupling constant. The generalized 
LY equation becomes 

8a 2 V 2 - a 2 R<j) + g~ l { n % t t: t J + 7 2 /3j#' )0- ? 

- r V W = 0. (21) 

This is just as well-behaved as the original LY equation: 
it always possesses a positive unique solution. This con- 
formal factor maps the free data onto a solution both of 
the generalized Hamiltonian constraint, Eq. ([20)) . and the 
momentum constraint. 

We can do slightly better than this. Let us assume 
that we would like to have a more general Hamiltonian, 
and add an R 2 term to the potential. Thus, we would 
like initial data which satisfies 

-a 2 gR-pgR 2 + g lk g 3l ii l ^ kl + 1 2 p)Pl-p 2 iT 2 = 0, (22) 

where p is another parameter. The analogue of the LY 
equation will be 

8a 2 V 2 - a 2 R4> + g- x + 7 2 /3^)<T 7 - Ptf^ 

= p^^ij 2 ^- R<p) 2 . (23) 

This, as mentioned earlier, is a deeply unpleasant equa- 
tion. However, let us linearize it about p = 0. First, let 
us solve Eq. (|2"Tt and use that conformal factor to map the 
free data to a set that satisfies Eq. (f22j) . This means that 
we are varying about <fro = 1. Then, differentiate Eq. ([23|) 
with respect to p, and set p = and tj> = 1. Denoting 
d(f)/dp — h, we get 

8a 2 W 2 h~a 2 Rh~7g-\Tr^ T ^ T +j 2 ^fi 3 i )h~5l3 2 p 2 h = R 2 . 

(24) 

Note that the h terms that should appear on the right 
hand side vanish because they are all multiplied by p. 
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Multiplying Eq. ([20l) by h/g and subtracting from Eq. (j24l 
(to eliminate the term linear in R) we get 



8a 2 V 2 h- (8g 



- 1 

TT, 



TT"ij 



A(3 2 p 2 )h = R 2 . 

(25) 

It is obvious that the coefficient of the undifferentiated h 
in Eq. (|2"5")) is negative. This means that the homogeneous 
equation has no kernel. In turn, this means that the 
inhomogeneous equation, Eq. (1231) , has a unique solution. 
This is the Fredholm alternative (see, for example 
Now we can use the implicit function theorem (again, see 
Hi) to guarantee that the non- linear equation, Eq. ([23| . 
has a solution for a range of p's in a neighborhood of 
zero. Unfortunately, this technique gives us no estimate 
as to the size of this neighborhood. One can immediately 
see that this technique works for any choice of metric 
potential. 

We construct CMC initial data, not because it plays 
any fundamental role in this theory, but simply out of 
convenience. The slices generated by the Hamiltonian, 
Eq. (|14p . will not stay CMC. We can use the same tech- 
nique, i.e., the Fredholm alternative -I- the implicit func- 
tion theorem, to relax the condition that p is a constant. 
We can replace it by assuming p = po + Opi , where po is 
a nonzero constant, p\ is a function of position, and 9 is 
a parameter. Now the conformal method gives us a sys- 
tem of 4 coupled nonlinear equations 12]. If we linearize 
about 9 = 0, we hnd that the equations decouple and we 
get a system with the necessary existence and uniqueness 
properties. This means we can construct non-CMC ini- 
tial data, but it is not clear how large the deviation from 
CMC we can allow. 

We can add a whole range of matter fields to the sys- 
tem. We can add a cosmological constant, a massive 
or massless scalar field, a Maxwell or a Yang-Mills field, 
dust, or even neutrinos The reader should note that 
the coupling constant of the Cotton- York density term in 
the generalized Hamiltonian, r y 2 (3j(3f, is dimensionless in 
natural units (in Eq.(l) and elsewhere in this article we 
use units 2k=1). In 3-covariant modifications of general 
relativity this term is essential to the perturbative power- 
counting renormalizability of the quantum theory [2], [H| ■ 

At early (intrinsic) times, the R term in Eq.(17) is 
suppressed by the e lnfl factor, while the Cotton- York ten- 
sor density-squared term is conformally invariant and in- 
dependent of g; and at late times the theory becomes 
more and more like GR. Thus, near the Big Bang, with 
ln<? — > — oo, the Cotton- York density term should domi- 
nate. In particular, we expect that the BKL model [lBT ]. 
where the time variation of the gravitational field is ex- 
pected to dominate over the spatial variations, will no 
longer hold in the presence of non-trivial amount of the 
Cotton- York tensor. 



A key advantage the intrinsic time formalism has over 
the extrinsic time formalism is that the Hamiltonian con- 
straint can be easily put in the form of an algebraic equa- 
tion for tt. If we wanted to use extrinsic time, we would 
need to solve the Hamiltonian constraint for ^fg. There- 
fore the reduced Hamiltonian in the extrinsic time gauge 
is a non-local object. Another sign of non-locality is that 
the lapse function for a CMC foliation is determined by 



V 2 N 



K{t). 



(26) 



This is obviously a nice elliptic equation, but clearly non- 
local as distinct from Eq. (ll3[) . the equation for the intrin- 
sic time lapse. 

It is clear that using an intrinsic time gives a very clean 
Hamiltonian structure in classical gravity. These good 
properties carry over when we implement the canonical 
quantization program Q. 
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